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Abstract —In the present analysis, we study the two-dimensional 
unsteady MHD free convection laminar heat and mass transfer 
flow past a semi-infinite vertical porous plate embedded in a 
porous medium in a slip flow regime in presence of thermal 
radiation, heat source/sink and Soret effects. A magnetic field of 
uniform strength is assumed to be applied transversely to the 
direction of the main flow. Perturbation technique is applied to 
transform the non-linear coupled governing partial differential 
equations in dimensionless form into a system of ordinary 
differential equations. The equations are solved analytically and 
the solutions for the velocity, temperature and concentration 
fields are obtained. The effects of various flow parameters on 
velocity, temperature and concentration fields are presented 
graphically. For different values of the flow parameters involved 
in the problem, the numerical calculations for the skin-friction 
coefficient, Nusselt number and Sherwood number at the plate 
are performed in tabulated form. 

Index Terms — MHD, skin-friction coefficient Introduction, 
Magnetohydrodynamic. 

I. INTRODUCTION 

Natural convection flow over vertical surfaces 
immersed in porous media has paramount importance because 
of its potential applications in soil physics, geohydrology, and 
filtration of solids from liquids, chemical engineering and 
biological systems. Ostrach [1], the initiator of the study of 
convection flow, made a technical note on the similarity 
solution of transient free convection flow past a semi-infinite 
vertical plate by an integral method. Study of fluid flow in 
porous medium is based upon the empirically determined 
Darcy’s law. Such flows are considered to be useful in 
diminishing the free convection, which would otherwise occur 
intensely on a vertical heated surface. Study of flow problems 
through porous medium is heavily based on Darcy’s 
experimental law [2]. Wooding [3] and Brink- man [4,5] have 
modified Darcy’s law, which are used by many authors on 
study of convective flow in porous media. 

Many transport processes exist in nature and in industrial 
applications in which the simultaneous heat and mass transfer 
occur as a result of combined buoyancy effects of thermal 
diffusion and diffusion of chemical species. 
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A few representative fields of interest in which combined 
heat and mass transfer plays an important role are designing 
of chemical processing equipment, formation and dispersion 
of fog, distribution of temperature and moisture over 
agricultural fields and groves of fruit trees, crop damage due 
to freezing, and environmental pollution. In this context, 
Soundalgekar [ 7] extended his own problem of Soundalgekar 
[6 ] to mass transfer effects. Callahan and Marner [8] 
considered the transient free convection flow past a semi¬ 
infinite vertical plate with mass transfer. Unsteady free 
convective flow on taking into account the mass transfer 
phenomenon past an infinite vertical plate was studied by 
Soundalgekar and Wavre [9], 

Magnetohydrodynamic flows have many 
applications in solar physics, cosmic fluid dynamics, 
geophysics and in the motion of earth’s core as well as in 
chemical engineering and electronics. Verma and Srivastava 
[10] proposed the effect of magnetic field on unsteady blood 
flow through a narrow tube. Huges and Young [11] gave an 
excellent summary of applications. MHD free convective 
flow of an electrically conducting fluid between two heated 
parallel plates in the presence of induced magnetic field by 
using analytical solution was carried out by Sharma [12]. 
Chakrabarti and Gupta [13 ] explained the heat transfer effect 
on hydromagnetic flow over a stretching sheet. Convective 
heat transfer effect on MHD flow past a continuously moving 
plate embedded in a non-Darcian porous medium has been 
proposed by Abo-Eldahab and El-Gendy [14 ]. 

Radiation is a process of heat transfer through 
electromagnetic waves. Radiative convective flows are 
encountered in countless industrial and environmental 
process. For example, heating and cooling chambers, fossil 
fuel combustion energy processes, evaporation from large 
open water reservoirs, astrophysical flows etc. Radiative heat 
and mass transfer play an important role in space related 
technology. The effect of radiation on various convective 
flows under different conditions have been studied by many 
researchers including Hussain and Thakar [15], Ahmed and 
Sarmah [16], Rajesh and Varma [17], Pal and Mondal [18], 
Samad and Rahman [19], Karthikeyan et al. [20], Das et al. 
[21] and Pal ef a/. [22], 

The study of heat source/sink effects on heat transfer 
is very important because its effects are crucial in controlling 
the heat transfer. Radiation effects on unsteady MHD flow 
through a porous medium with variable temperature in 
presence of heat source/sink is studied by Vijaya Kumar et al. 
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[23], Vijaya sekhar and Viswanadh reddy [24 ] have obtained 
the analytical solution for the effects of heat sink and 
chemical reaction on MHD free convective oscillatory flow 
past a porous plate with viscous dissipation. Gireesh Kumar 
and Satyanarayana [25] studied the heat and mass transfer 
effects on unsteady MHD free convective waiter’s memory 
flow with constant suction in presence of heat sink. 
Choudhury and Paban Dhar [26] investigated the effects of 
MHD Visco-elastic fluid past a moving plate with double 
diffusive convection in presence of heat 
generation/absorption. 

At the macroscopic level, it is accepted that the 
boundary condition for a viscous fluid at a solid wall is one of 
“no slip”. While no-slip boundary condition has been proven 
experimentally to be accurate for a number of macroscopic 
flows, it remains an assumption that is based on physical 
principles. In fact nearly two hundred years ago, Navier’s [27] 
proposes a general boundary condition that incorporates the 
possibility of fluid slip at a solid boundary. Navier’s proposed 
condition assumes that the fluid slip velocity at a solid surface 
is proportional to the shear stress at the surface. The 
mathematical form of the Navier’s proposed condition on slip 


f dv^ 

velocity as emphasized by Goldstein [28] is V = y — 

V d yj 

where y being the slip coefficient, v the slip velocity and y 


the normal coordinate. The fluid slippage phenomena at solid 
boundaries appear in many applications such as in micro¬ 
channels or nano channels and in applications where a thin 
film of light oil is attached to the moving plates or when the 
surface is coated with special coating such as thick monolayer 
of hydrophobic octadecyl- trichorosilane. Due to practical 
applications of the fluid slippage phenomenon at solid 
boundaries, several scholar have carried out their research 
work in that literature, the names of whom (Yu & Ahem [29], 
Waltanebe et al. [30], Jain and Sharma [31], Khaled and Vafai 
[32], Saxena and Dubey [33 ] and Poonia and Chaudhary 
[34]) are worth meaning. 


Due to the importance of Soret (thermal-diffusion) 
and Dufour (diffusion-thermal) effects for the fluids with very 
light molecular weight as well as medium molecular weight. 
The Soret effect arises when the mass flux contains a term 
that depends on the temperature gradient. Bhavana et al [35] 
proposed the Soret effect on unsteady MHD free convective 
flow over a vertical plate in presence of the heat source. 
Anand Rao et al. [36] explained the Soret and Radiation 
effects on unsteady MHD free convective flow past a vertical 
porous plate. Soret effect on MHD flow of heat and mass 
transfer over a vertical stretching plate in a porous medium in 
presences of heat source was proposed by Mohammad Ali and 
Mohammad Shah Alam [37], Dufour and Soret effects on 
unsteady MHD free convective flow past a vertical porous 
plate embedded in a porous medium with mass transfer was 
studied by Alam et al [38 ]. Mohammed Ibrahim [39] studied 
the Soret and Dufour effects on unsteady MHD convective 
heat and mass transfer flow past an infinite vertical plate 
embedded in a porous medium in presence of radiation. 

The present work is concerned with the effect of 
mass transfer and thermal radiation on magnetohydrodynamic 


convection flow of an unsteady viscous incompressible 
electrically conducting fluid past a semi-infinite vertical 
porous plate immersed in a porous medium in slip flow 
regime in presence of heat sink and Soret effects. The 
classical model for radiation effect introduced by Cogley et al. 
[40] is used. Perturbation technique is applied to convert the 
governing non-linear partial differential equations into a 
system of ordinary differential equations, which are solved 
analytically. 


II. MATHEMATICAL ANALYSIS 


x 



Figure 1. Physical model of the problem. 

We consider a two-dimensional unsteady flow of a 
laminar, incompressible, electrical conducting and heat 
absorbing fluid past a semi-infinite vertical porous plate 
embedded in a uniform porous medium. We introduce the 

coordinate system (x, y, Z ) with X axis is chosen along the 

plate, Y axis perpendicular to it and directed in the fluid region 
and Z axis is along the width of the plate as shown in the 

Figure 1. A uniform magnetic field of strength B () in the 

presence of radiation is imposed transversely in the direction 
of Y axis. The induced magnetic field is neglected under the 
assumption that the magnetic Reynolds number is small. It is 
assumed that there is no applied voltage which implies the 
absence of any electrical field. The radiative heat flux in the X 
direction is considered negligible in comparison to that in Y 
direction. The governing equations for this study are based on 
the conservation of mass, linear momentum, energy and 
species concentration. Taking in to consideration the 
assumptions made above, these equations in Cartesian frame 
of reference are given by equation of 

Continuity: 



( 1 ) 


Momentum equation: 
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du _du 
—+v — 
dt dt 


1 dp 
pdx 


S/>(f-?,)+s/(C-C,)+v 


, 2 - 


du 

w 


cB-_ v 

- -U -r 

D K 


( 2 ) 


where / = K Zw —d A, K /w is the absorption 

J " QJ 

coefficient at the wall and e hA is the Planck’s function. 


Energy equation: 


dT_ _1 dg r Qo(T -Tj) 

8T V 8T pC p dy 2 pC p dy pC p 

(3) 


Species equation: 


8C _ 8C 
— + v — 
dt dt 


= D 


M 


8 2 C _ 8 2 T 

+ ZX 




w 


(4) 


Under the assumption, the appropriate boundary conditions 
for velocity involving slip flow, temperature and 
concentration fields are given by 

_ PjJT _ _ _ _ * * 

« = =h—, T = T w + s(T„-T,)e"' , 

C = C„ + £(C„-C„)e" V a, j=0 (6) 

i7^£7.=t/„(l + ffe" v ), C ^ C as 

y —> oo (7) 


where x , y and t are the dimensional distances along and 
perpendicular to the plate and dimensional time, respectively. 
U and V are the components of the dimensional velocities 
along X and y respectively, p is the density of the medium, 
g is the acceleration due to gravity, V is the kinematic 
viscosity, (7 is the fluid electrical conductivity, B 0 is the 

magnetic induction, K is the permeability of the porous 
medium, [3 is the coefficient of thermal expansion, [3 is the 
coefficient of mass expansion, T is the dimensional 
temperature of the fluid near the plate, T w is the dimensional 

free stream temperature, C is the dimensional concentration 

of the fluid near the plate, C x is the dimensional free stream 

* 

concentration, k is the thermal conductivity of the fluid, q r is 
the radiative heat flux, the term Q {) (T — T x ) is assumed to 
be amount of heat generated or absorbed per unit volume, Q 0 
is constant, which may take on either positive or negative 
values. When plate temperature T exceeds the free stream 
temperature T x , the source term Q {] > 0 and heat sink when 
< 0, // is the fluid viscosity, C is the specific heat at 
constant pressure, D M is the coefficient of chemical 
molecular diffusivity and D T is the coefficient of thermal 
diffusivity. 

Cogley et al. [40 ] showed that, in the optically thin limit for a 
non-gray gas near equilibrium, the radiative heat flux is 
represented by the following form: 

^L = a(t-T\ i* 

dy 1 ’ 

(5) 


where T and C w are the dimensional temperature and 

species concentration at the wall respectively and ll is the 
characteristic dimension of the flow fluid. The suction 
velocity normal to the plate is a function of time only, it can 
be taken in the exponential form as 


v =-V 0 (l + sAe nY J 


( 8 ) 


where A is a real positive constant, £ and £ A are small 
quantities less than unity and V 0 is a scale of suction velocity 
which is a non-zero positive constant. 

In the free stream, from equation (2) we get 


1 8p 

p dx 


dU , aB: 


dt 


. + ^LU + ^U 

oo T r* a 

P K 


(9) 


Now we introduce the dimensionless variables as follows 


U V yV 0 

u = —, v = — , y = —- 
U n K v 


U = 


u n 


t = 


tyi 

V 


T-T 
0 = = 
T-T 


</> = 


c-c„ 

c.-c 
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n v 
V 2 

v n 


KV 
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v 


Pr = 


P C B 
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( 10 ) 
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Gr = 


Q = 


vfig(T w -TJ 


U V 

u o y o 


Gm = 


vfg(C w -CJ 


U V 

^ o v o 


Qq v 


However, these equations can be reduced to a set of ordinary 
differential equations, which can be solved analytically. These 
can be done by representing the velocity, temperature and 
concentration of the fluid in the neighborhood of the plate as 


R = 


pc p v 0 2 

4 vf 

PCX 


V h 
h = _2_ 


V 


So = 

V 


r ^'\V ^co 

C -C 

V'-'w '“'oo 


Sc = 


v 


D 


M 


u = u Q (y) + se nt u x (y) + OX ) 
9=9 0 (y) + se m 0 l (y) + O(s 2 ) 
^ = (4 (y) + ee m fa(y) + 0(s 2 ) 


(16) 


(17) 


(18) 


where Pr is the Prandtl number, M is the magnetic field 
parameter, Gr is the Grashof number for heat transfer, Gm is Substituting (16) to (18) in equations (11) to (13) and equating 
the Grashof number for mass transfer, Q is the heat sink the harmonic and non harmonic terms and neglecting the 

parameter, a is the permeability parameter, 6 is the non coe ffi c ient of 0(£ 2 ) we get the following pairs of equations 

dimensional temperature, (j) is the non dimensional 
concentration, R is the radiation parameter, h is the rarefaction 
parameter. So is the Soret number and Sc is the Schmidt 
number. 


for (u 0 ,9 0 ,<f) 0 ) and (u ] , 67, . 

u" + u' 0 — Nu 0 = -N - Gr6 {) — Gm(/) l: 


In view of equations (8) to (10) the governing equations (2), 
(3) and (4) reduce the following non-dimensional form: 

* (1 

v ’ dy dt dy- 


(19) 


u"+u[ — (N + n)u l = —Au' 0 — GrO x — Gmfy — (N + n ) 

( 20 ) 


dt 


de_ 

dt 

d(j) 

dt 


(l + sAe n 


\d0 d-d 
^ dy Pr dy 2 


( 11 ) 


-R9-Q6 


d 9 


(l + £ Ae nt )^- = — ^4- +So 
' ’ dy Sc dy " dy 2 


( 12 ) 


(13) 


9” + Pr 9' } - Pr(R + Q)9 0 = 0 

9"+ Pr O' - Pr (R + Q + n )9 { = -A Pr 9' 

(f>"+Sc(/); ) =-scSoO" 

(/>"+ Sc(f>[ - nSc(J\ = -ASc(/)' q - ScSoO" 


( 21 ) 


( 22 ) 


(23) 


(24) 


where N = M -\ - 

K 

The boundary conditions (6) and (7) in the dimensionless 
form can be written as 

u = u u =h —, 9 = 1 + se n ' , (f> = \ + se nt at y = 0 


dy 


(14) 


where the primes denote the differentiation with respect to y. 
The corresponding boundary conditions can be written as 

u Q = hu' 0 , u x = hu[ , 9 0 = 1, 9 X = 1, (f) Q = 1, (f> x = I at y 
= o (25) 

My = 1, Mj = 1, Oq —^ 0, 9 X —^ 0, <f>Q —^ 0, cj) x —^ 0 as 
y —>• qo (26) 


U —> U x — 1 + se'" , 9 0, ^ > 0 as y — >• oo The solutions of equations (19) to (24) which satisfy the 

(15) 


III. SOLUTION OF THE PROBLEM 

Equations (11) to (13) are coupled non-linear partial 
differential equations and these can be solved in closed form. 


boundary conditions (25) and (26) are given by 

u 0 = 1 + A l0 e-’^ y - A n e nhy - \e ,rhy + A^e^ y 


( 21 ) 
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-'» 2 .v Knowing the temperature field, it is interesting to study the 

non-dimensional form of the rate of heat transfer in terms of 
Nusselt number at the plate is given by: 


«, =1+- 4<f v + V miy -A 15 e 


-A l6 e~ m ' y - \ 7 e m * y - - A 20 e~ m ' } ’ 


(28) 


e 0 =e- ,n ' y 


6» =\e~" hy +\e 


-m 3 y 


— A 



(29) 

r m \y 

(30) 

-A 2 e-' n ' y 

(31) 

’- m3y +A 4 e- nby +A 5 e~ m ' y 

(32) 


In view of the above solutions, the velocity, temperature and 
concentration distributions in the boundary layer become 


N„ =- 


f d0 A 


so, 


o 

mj + ce" (-m 2 A () — in l A 1 )j 


0 +se nt 86 ' 




Knowing the concentration field, it is interesting to study the 
non-dimensional form of the rate of mass transfer in terms of 


S h = - 


r d^ 

V^A=o 


f d k +6 . e n, 


l=o 


-(—(1 + A 2 )m 3 + m l A 2 + se nt (—m 4 A (t — m 3 A 3 — m 2 A 4 — 


IV. RESULTS AND DISCUSSION 


m, 


u(y,t) =1+ -A,e~ v - +fe "(l++Ai^ - 

-V*+A/”' 3 ’-a/* - ■-■V~' v -As^-A/"” 2 ’’ - AoO 

6>(_y,0 = U %,! +£e'''(A^ %, +A^ mi)! ) 

fly,!)=(1+A)e"^ ^- Ae""'"+A e_v +A e " w +■ ¥' v ) 

It is now important to calculate the physical quantities of 
primary interest, which are the local wall shear stress, the 
local surface heat and mass flux. Given the velocity field in 
the boundary layer, we can now calculate the local wall shear 
stress (i.e., skin-friction) is given by 

* ( du ' 

T W = P ■= 

l^l=o 

and in dimensionless form, we obtain 


In order to get physical insight into the problem, we have 
carried out numerical calculations for non-dimensional 
velocity field, temperature field, concentration field, co¬ 
efficient of skin-friction C f - at the plate, the rate of heat 

transfer in terms of Nusselt number N and the rate of mass 
transfer in terms of Sherwood number S h by assigning 

specific values to the different values to the parameters 
involved in the problem, viz., Magnetic field parameter M, 
Grashof number for heat transfer Gr , Grashof number for 
mass transfer Gm, permeability parameter K, heat sink 
parameter Q, rarefaction parameter h, Prandtl number Pr , 
Radiation parameter R, Schmidt number Sc and Soret number 
So, time t. In the present study, the following default 
parametric values are adopted. Gr = 6.0, Gm = 4.0, M = 3.0, 
K = 1.0, n = 0.1, A = 1.0, t = 1.0, Pr = 0.71, R = 1.0, Q = 
1.0, Sc = 0.6, So = 1.0, h = 0.3 and £ = 0.2. All graphs 
therefore correspond to these unless specifically indicated on 
the appropriate graph. The numerical results are demonstrated 
through different graphs and table and their results are 
interpreted physically. 


C / = 


^ du^ 


pu 0 V 0 


d y J v=o 


= u\ 0) = 


8u t 

~dy 


° + se’ Su ' 


dy 


J v=0 


( y -m i A itj +m l A 1 +m 3 A i -m l \ j )+se'"(-m 6 A B -m 5 A n +m l A a +m } A D -m l A u +m 1 A ls + 

+m l A l6 + m 4 A l7 + m i A lg + m 3 A X9 + m l A 20 ) 


It is observed from Figure 2 that an increase in 
Grashof number for heat transfer Gr leads to a rise in the 
values of velocity u due to enhancement in buoyancy force. 

The plot of velocity profile for different values of 
Grashof number for mass transfer Gm is given in Figure 3. It 
is observed that velocity increase for the increasing values of 
Grashof number for mass transfer Gm. 

Figure 4 plots the velocity profiles against the span- 
wise coordinate y for different magnetic field parameter M. 
this illustrates that velocity decreases as the existence of 
magnetic field becomes stronger. This conclusion agrees with 
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the fact that magnetic field exerts retarding force on the free- 
convection flow. 

The changes in velocity profile due to different 
permeability of the porous medium are plotted in Figure 5. 
This figure shows that the velocity profiles increases rapidly 
as increases permeability of porous medium K. 

Figure 6. illustrates the dimensionless velocity u for 
different values of the Prandtl number Pr. The analytical 
results show that the effect of increasing values of Prandtl 
number number results in a decreasing velocity. 

The effect of the Radiation parameter R on the 
dimensionless velocity u is shown in Figure. 7 Figure7 shows 
that velocity component decreases with an increase in the 
radiation parameter R. 

Figure 8 depict the dimensionless velocity 
component u profiles for different values of heat sink 
parameter Q. It is noticed that an increase in the heat sink 
parameter Q results in decrease in the dimensionless velocity 
component u within the boundary layer. 

The influence of the Schmidt number Sc on velocity 
profiles are plotted in Figure 9. The Schmidt number Sc 
embodies the ratio of the momentum to the mass diffusivity. It 
is noticed that 

as the Schmidt number Sc increases the velocity decreases. 
This causes the concentration buoyancy effects to decrease 
yielding a reduction in the fluid velocity. 

Figure 10 depicts the velocity profiles for different 
values of the Soret number So. The Soret number defines the 
effect of the temperature gradients inducing significant mass 
diffusion effects. It is noticed that an increase in the Soret 
number results in an increase in the velocity within the 
boundary layer. 

Figure 11 indicated that the velocity component u for 
different values of rarefaction parameter h. The velocity it 
increases as the rarefaction parameter h is increased indicating 
the fact that slips at the surface accelerates the fluid motion. 

The horizontal velocity profile in the boundary layer 
for different values of time t is depicted in Figure 12. It is 
observed that the horizontal velocity slowly attains the peak 
value close to the porous boundary and then it decreases till it 
reaches the minimum value at the end of the boundary layer 


for all the values of time. It is noticed that the velocity 
increases as the time increases. 

From Figure 13. it is observed that an increase in the 
Prandtl number results a decrease of the thermal boundary 
layer thickness and in general lower average temperature 
within the boundary layer. The reason is that smaller values of 
Pr are equivalent to increasing the thermal conductivities, and 
therefore heat is able to diffuse away from the heated plate 
more rapidly than for higher values of Pr. Hence in the case 
of smaller Prandtl numbers as the boundary layer is thicker 
and the rate of heat transfer is reduced. 

Figure 14 shows the variation of temperature profiles 
with respect to the radiation parameter R. from this figure, it is 
observed that as temperature increases for the increasing 
values of radiation parameter R. This result qualitatively 
agrees with expectations, since the effect of radiation is to 
decrease the rate of energy transport to the fluid, thereby 
decreasing the temperature of the fluid. 

The influence of the parameter heat sink on 
dimensionless temperature profiles 9 is plotted in Figure 15. 
It is noticed that dimensionless temperature decreases with an 
increase in heat sink parameter Q. 

Moreover Figures 16 displays that the temperature 
profiles 9 increases as time t is increased. 

The concentration profiles for different values of 
Schmidt number Sc are plotted in Figure 17. For Schmidt 
number Sc, the value 0.6 corresponds to water vapor and 
represents a diffusing chemical species of most common 
interest in air. The analytical results show that the effect of 
increasing Schmidt number in air and water results in a 
decreasing concentration distribution across the boundary 
layer. The results also indicate that the effect on Schmidt 
number on the velocity and concentration is more in water as 
compared to that in air. 

For various values of Soret parameter So, the 
concentration profiles is plotted in Figure 18. Clearly as So 
increases, the dimensionless concentration increases. 

Concentration profiles (j) for different values of £ 
and time t are shown in Figures 19. From this figure, it is 
observed that dimensionless concentration profiles increases 
as time t is increased. 
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Figure 2. Velocity u versus y, under the effect of Gr 

Figure 3. Velocity u versus y, under the effect of 
Gm 
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Figure 4. Velocity u versus y, under the effect of M 


Figure 5. Velocity u versus y, under the effect of 


K 




y 


Figure 6. Velocity u versus y, under the effect of Pr Figure 7. Velocity u versus y, under the effect of 
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Figure 10. Velocity n versus y, under the effect of So 


Figure 11. Velocity u versus y, under the effect 
of h 
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Figure 12. Velocity u versus y, under the effect of t 
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Figure 13. Temperature 9 versus y, under the 
effect of Pr 
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Figure 14. Temperature 9 versus y, under the effect 
of R 



Figure 15. Temperature 9 versus y, under the 
effect of Q 




Figure 16. Temperature 9 versus y, under the effect 
of t 


Figure 17. Concentration (j) versus y, under the 
effect of Sc 
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Figure 18. Concentration (j) versus y, under the effect 
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Figure 19. Concentration (j) versus y, under the effect of t 


of So 
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Table 1. Skin friction coefficient, Nusselt number and Sherwood number for various values of Gr, Gm, M, K, h with Pr = 0.71, 
R = 1.0, Q = 1.0, Sc = 0.6, So = 1.0. £ = 0.2, n = 0.1,A = 1.0, t = 1.0. 


Gr 

Gm 

M 

K 

h 

Cf 

V, 

s„ 

6.0 

4.0 

3.0 

1.0 

0.3 

4.5564 

2.0577 

0.1785 

2.0 

4.0 

3.0 

1.0 

0.3 

3.7052 

2.0577 

0.1785 

4.0 

4.0 

3.0 

1.0 

0.3 

4.1308 

2.0577 

0.1785 

6.0 

1.0 

3.0 

1.0 

0.3 

3.4571 

2.0577 

0.1785 

6.0 

2.0 

3.0 

1.0 

0.3 

3.8235 

2.0577 

0.1785 

6.0 

4.0 

1.0 

1.0 

0.3 

5.5261 

2.0577 

0.1785 

6.0 

4.0 

2.0 

1.0 

0.3 

4.9226 

2.0577 

0.1785 

6.0 

4.0 

3.0 

2.0 

0.3 

4.7198 

2.0577 

0.1785 

6.0 

4.0 

3.0 

3.0 

0.3 

4.7823 

2.0577 

0.1785 

6.0 

4.0 

3.0 

1.0 

0.5 

3.5121 

2.0577 

0.1785 

6.0 

4.0 

3.0 

1.0 

0.7 

2.8571 

2.0577 

0.1785 


Table 2. Skin friction coefficient, Nusselt number and Sherwood number for various values of Pr, R, Q, Sc, So with Gr = 6.0, 
Gm = 4.0, M = 3.0, K= 1.0, h = 0.3, £ =0.2, ft = 0.1, A = 1.0, t= 1.0. 


Pr 

R 

Q 

Sc 

So 

C / 

V, 

s h 

0.71 

1.0 

1.0 

0.6 

1.0 

4.5564 

2.0577 

0.1785 

1.0 

1.0 

1.0 

0.6 

1.0 

4.4364 

2.5951 

- 0.0256 

3.0 

1.0 

1.0 

0.6 

1.0 

4.0611 

5.8503 

- 1.1499 

0.71 

2.0 

1.0 

0.6 

1.0 

4.4873 

2.3689 

0.0337 

0.71 

3.0 

1.0 

0.6 

1.0 

4.4350 

2.6350 

- 0.0869 

0.71 

1.0 

2.0 

0.6 

1.0 

4.4873 

2.3689 

0.0337 

0.71 

1.0 

3.0 

0.6 

1.0 

4.4350 

2.6350 

- 0.0869 

0.71 

1.0 

1.0 

0.78 

1.0 

4.4478 

2.0577 

0.3835 

0.71 

1.0 

1.0 

0.94 

1.0 

4.3351 

2.0577 

0.6758 

0.71 

1.0 

1.0 

0.6 

2.0 

4.7931 

2.0577 

-0.5107 

0.71 

1.0 

1.0 

0.6 

3.0 

5.0298 

2.0577 

- 1.1998 
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